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Abstract 

We present a brief discussion on the nonlinear Schrodinger equa- 
tion for modeling the propagation of the deep-water wavetrains and a 
discussion on its doubly-localized breather solutions that can be con- 
nected to the sudden formation of extreme waves, also known as rogue 
waves or freak waves. 
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1 Introduction 



Nonlinear phenomena are much studied in various area of science pQ [6]. 
Many of these phenomena are modeled by nonlinear partial differential equa- 
tions or systems of such equations, which increased the interest in the meth- 
ods of obtaining exact and approximate solutions of such equations [7)-[TT]. 
In this brief report we shall be interested in extreme deep-water phenom- 
ena, also referred to as freak or rogue waves [T21 [13] and its description by 
breather solutions of the nonlinear Schrodinger equation. Rogue water waves 
are extreme high sea waves that can cause severe damage on commercial and 
other ships or on oil platforms. Such waves are deep-water waves which prob- 
ably can be described by breather solutions of equations that belong to the 
family of the nonlinear Schrodinger equation. The latter will determine the 
scope of our report as follows: (I) short notices on linear and nonlinear sur- 
face gravity water water; (II) Heuristic derivation the nonlinear Schrodinger 
equation (NSE) for deep-water waves; (III) discussion on one extension of 
this equation (known as Dysthe equation); and (IV) a short discussion on 
breather solutions of the NSE obtained by Peregrine as well as Akhmediev 
and coworkers , which could be considered as appropriate models for the 
freak waves. 



2 Model equations for water waves: incom- 
pressible inviscid approximation 

We start from the differential form the basic equation of fluid mechanics are 
as follows [T7j 

: mass conservation, 

1 77* 
— \7p + f H V 2 v : momentum conservation. 

P * P 

In EqsQ p is the density of the fluid; v is the fluid velocity; p is the pressure; 
/ summarizes the body forces acting on the fluid; rj* is coefficient of viscosity 
called kinematic viscosity. 

In addition boundary conditions should be imposed. We assume that the 
depth of the fluid is h and it is bounded from below by a hard horizontal 
bed. The upper fluid surface is assumed to be free. The unperturbed free 
upper surface is at z = 0. When the upper surface is perturbed, there 
is vertical displacement rj(x,y,t) of each point of the surface. Then, the 



(1) f t +(v-V)v 
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boundary condition on the upper fluid surface is at z = r){x,y,t). On the 
lower solid surface the normal component of velocity has to vanish, i.e. no 
flux is permitted at the bottom. That is, v z = at z = — h. 

We shall discuss the case of constant fluid density p = const (incompress- 
ible fluid approximation). In addition, the waves will have wavelength much 
longer than approximately 1.8 cm. For this case the viscosity effects are 
negligible, i.e., we shall consider gravity waves (and not capillary waves). 



2.1 Irrotational approximation and model of small am- 
plitude surface gravity waves 

The quantity uj = V x v is called vorticity of the flow and when uj = the 
flow is called irrotational. In this case the velocity of the flow is a potential 
field: v = V0, where (ft is the velocity potential (the central quantity we 
shall discuss below). In addition, for gravity waves / = —ge Z) where g 
is the acceleration of gravity and e z = (0,0,1). For small water surface 
amplitudes and in the irrotational approximation, the model equations ([!]) 
and the boundary conditions are reduced as follows 



0; —h < z < r)(x,y,t) —> from mass conservation, 

2 / o / \ 2 / o , \ 2> 



dt 

d(f) 
dz 



dz 




dz 



rig 



; at z = rj(x, y ) t) — )► from 



\dy J 

momentum conservation, 
dr\ dcj) drj d(j) dr] d(j) 

^r^r + ^r^r + ^r^r; at z = v{ x ,y,t) -> from b. con 

ox ox oy oy oz oz 
the top surface, 

0; at z = —h —> from b.c. on the bottom surface. 



(2) 



Thus, the model equation becomes linear (the Laplace equation) but the 
boundary conditions are nonlinear. 



2.2 Shallow- water waves and deep-water waves 

For small amplitudes (but long wavelengths) water waves, the nonlinear re- 
lationships from the system ^ can be linearized. If the mean surface dis- 
placement and the mean velocity potential are small with respect to the wave- 
length and to wave period scales, then, the nonlinear terms in the boundary 
conditions can be neglected. After a Taylor series expansion of the small 
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quantity 77 (and keeping only the first term of the expansion) the top bound- 
ary condition can be written as condition on z = 0. Thus, we obtain the 
following simplified (and linear) problem 



V 2 S = 0; 



(3) 



dt 



dz 



-g- 



dz 
0; at z 



h < z < 0, 

6 

at z — 0, 



-h. 



The next approximation is that we assume that the waves propagate in the 
x-direction and are uniform in the ^/-direction. Thus, the problem becomes 
one-dimensional and one searches for a traveling wave solution with wave 
frequency uj and wavenumber k: 



(4) 



0(x, t) = A{x, z) s'm(kx — ujt). 



The substitution of Eq. ^ in ^ leads to the following solutions for the 
velocity potential (j) and for the surface displacement r\ 



Tj = Acos(kx — u)t)\ (j) 



(5) 



. cosh(fc(z + h)) . /7 N 

ujA — — — — - — smlkx — out)] 

ksmh(kh) 

^ exp(-kh) sinh(/c/i), 



A 



where a is a constant of integration, A, k and uj denote the wave ampli- 
tude, the wavenumber and the wave frequency, respectively. The dispersion 
relation for the small amplitude surface water waves as well as their phase 
velocity v and group velocity v g are as follows 



1 + 



2kh 



2sinh(2fc/i) 



(6) uj 2 = gffctanh(fc/i); v = y — tanh(/c/i); v g 
The relationship R = de f th - r = - has two limit cases: R « 1 (shallow- 

^ wavelength A V 

water waves) and R » 1 (deep-water waves). For the case of shallow- water 
waves the dispersion relation ^ can be approximated by 

k 2 h 2 



UJ 



kJ gh 



1 
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+ ... 



c 



gh. 



kn 



For very long shallow water waves uj fcc ; v = ^ ~ cq; v g — || ~ #io 
deep-water waves the approximation for the dispersion relation is 
and the phase and group velocity are v 



For 



sj~gjk\ v g = yjg/(2k), respectively. 
That is, the group velocity is smaller and half the phase velocity. 
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3 Deep-water waves. The nonlinear Schrodinger 
equation. The Dysthe equation 

A weakly nonlinear model for shallow-water waves can be described by the 
Korteweg-de Vries equation. We shall be interested in the deep-water case, 
which can be modeled by the nonlinear Schrodinger equation. 

3.1 Derivation of the nonlinear Schrodinger equation 
by applying a Taylor series expansion to the dis- 
persion relation for deep-water waves 

A weakly nonlinear approximation to the nonlinear deep-water wave prob- 
lem are the Stokes waves. They however are unstable against modulation 
perturbations. The velocity of the Stokes waves to second-order in steepness 
is 

(7) »Wf< 1 + ^>- 

where a denotes now the wave amplitude. From Eq. ([7| one obtains easily the 

dispersion relation uj = \J gk{l + ^-). Let us consider a slowly modulated 
Stokes wave wavetrain 

(8) rj = Re[A(X,T)ex-p(i(u t - fax))], 

where ujq and fa are the frequency and wave number of carrier Stokes wave 
and A(X, T) is the modulation amplitude of the wavetrain. In addition, 
X = ex and T = et ( e « 1) are the slowly varying space and time 
variables, respectively. Physically, e := A fa is the steepness of the wave 
and is assumed to be small. Let us now perform a Taylor series expansion 
around the wavenumber fa and the amplitude A = A(0, 0)... The dispersion 
relation of the carrier Stokes wave is 



(9) u = ^gk(l + k 2 | A | 2 ), 

where | A | is the amplitude of the Stokes wave (and the amplitude of the 
envelope). The Taylor series expansion about the wavenumber fa of the 
carrier wave and about the envelope amplitude A = A = as follows [T7j 

(10) u = Uo + —{k - ko) + - w (k - k ) 2 + jy^jd A\ 2 - | Ao | 2 ). 

Let Q = uj — ujq and K = k — fa. In addition (accounting also for Eq.Q) 
|H \k=k = v g = ^ (note that the group velocity of the envelope is twice 
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smaller than the phase velocity of the carrier wave); \k=k — 2P = — 



Q = gpp Uo=o= ^ k 2 . Then, from Eq.(|lOf 
(11) tt = v g K + PK 2 + Q \A\ 2 . 

The Fourier and the inverse Fourier transform of the envelope function are 

/oo 
dXdT A(X, T) exp[i(QT - KX)], 
-oo 



A(X,T)=T- 1 [A(K,n)] 
(12) 



(I) 



/OO 
dKdtl A(K, Q) exp[-i(QT - KX)}. 
-oo 



From Eqs. (|12|) 
(13) 



-mF- l [A(K,Q)]. 



Q and K are of order e. Then from Eq.(13) we can write 

. d 



(14) 



K 



o • 9 

- te dx ; Q = te W 



The substitution of the relationships from Eq.(14) in Eq.(ll) and application 
of the resulting operator equation to the envelope amplitude A leads to the 
nonlinear Schrodinger equation for the evolution of the amplitude of the 
envelope of the wavetrain (e is incorporated in T and X by appropriate 
rescaling) . 



(15) 



'OA u dA^ 
df + 2k~ dX 



coo d 2 A _ 1 
8fcg dX 2 ~ 2 



u kl \A\ 2 A = 0. 



UJQ_rp 
2k ± 



Eq. (15) can be rescalled as follows: r = — fpT; £ = X — v g T = X 

(coordinate is in a frame that moves with the group velocity of the wavetrain); 
q = \/2fcg A The rescalled form of the nonlinear Schrodinger equation is 



(16) 



.do d 2 q rt , l9 



3.2 More accurate deep-water wave envelope equation: 
The Dysthe equation 

The Dysthe equation [18J is an extension of the nonlinear Schrodinger equa- 
tion. This equation aims to solve the problem with the bandwidth limitation 
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of the NSE. The nonlinear Schordinger equation is valid for small steepness 
values, i.e. k A « 1 and when the bandwidth is narrow (Ak/k « 1, 
Ah is the modulation wavenumber). In order to obtain the Dysthe equation 
one starts with the model equations for the velocity potential (j)(x,y,z,t) 
and surface displacement rj(x,y,t) for an incompressible inviscid fluid with 
uniform depth h 

V 2 (j) = 0, for - h < z < 77, 

dh ^ ^ 7 d(f) 

(17) ^ = 0, atz = -/i. 

Next one assumes that fcA = 0(e); ^ = 0(e) 1 / 2 , (fc/i)" 1 = C^e 1 / 2 ) and 
performs the following expansion for the velocity potential <fi and surface 
displacement h 

= ^+1 (Aexp(i6 + kz) + A 2 exp(2(i6» + kz)) + . . . + c.c.) , 



h = h+^{B exp(i6» + kz) + B 2 exp(2(i6» + fcz)) + . . . + c.c.) , 



(18) 



where 9 = kx — cut. The drift 0, the set down h and the amplitudes 

A, t4 2 , • • • ? -B, £>2, • • • are functions of the slow modulation variables ex, ey, ei. 
Furthermore, a non-dimensionalization of the variables is performed as fol- 
lows: ut -> t\ k{x,y,z) -> (x,y,z)', /c 2 cj _1 (,4, . . . , A n) 4>) -^_(A, . . . , A n , 0); 
k(B, . . . , B n , 77) — >> (£?, . . . , B n , 77). Note that the functions 0; 77; A, . . . , ^4 n ; 

B, ...,5 n are functions of the variables e 1 / 2 (x,z,t). Proceeding with the 
expansions up to order 0(e 3 5 ), one obtains the Dysthe equation 

<M 1 <M z d 2 A _id 2 A _ 1 <9 2 A 3 d 3 A H d A A 
~dt + Ydi + 8 d^ 2 " ~ 4 ~dyi ~ 16 + 8 5^^ ~ 128 5^ + 
15i d 4 A 37 d 4 A 7 2 7 <9 5 ,4 35 <9 5 ,4 21 <9 5 ,4 
32 dx 2 dy 2 ~ 32~d^ + 2 ' ' + 256 ~ 64 <9x 3 <97/ 2 + 64 dxd?/ 4 + 

(19) 2^1 ar + < = °- 

for the corresponding boundary conditions 

V 2 ^ = 0, for - h< z < 0, 
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(20) 



— = 0. at z — —h. 



4 Peregrine and Akhmediev-Peregrine breathers 

Akhmediev and co-workers [HI [15] derived a family of space-periodic pulsat- 
ing solutions of the NSE, also called Akhemdiev breathers, that start from 
the plane wave solution at r = — oc and return again to the plane wave form 
at r = oo. Taking this period to infinity, we obtain the Peregrine breather 
[16], a solution which is doubly-localized, in space and time, and pulsates 
only once. These properties make the Peregrine breather an ideal model 
to describe oceanic rogue waves. In the same work [2], Akhmediev and 
co-workers derived a higher-order doubly-localized breathers, also referred 
to the Akhmediev-Peregrine breather. In fact, there is an infinite hierar- 
chy of doubly-localized Akhmediev-Peregrine solutions. Generally, the j—th 
Akhmediev-Peregrine solution can be written in terms of polynomials: 

(21) qj(£, r) = q exp(2i | q | 2 r) 

where q is proportional to the carrier wave amplitude and r); Hj(£, r) 
and Dj(^t) are appropriate polynomials. The first-order rational solution 
(Peregrine breather) is given by 

(22) G x = 4; H, = 16 | q | 2 r; D 1 = 1 + 4 | q | 2 £ 2 + 16 | q | 4 r 2 

The second-order rational solution (Akhmediev-Peregrine breather) is given 
by 

G 2 =(k| 2 e 2 + 4ko| 4 r 2 + £) (|gor^ 2 + 20Uo| 4 T 4 + £)-^ 
H 2 = 2 | g | 2 r(4 \ q | 4 r 2 - 3 | q | 2 ^ 2 ) + 2 | q | 2 r ((2 | q | 2 e + 

4ko| 4 r 2 ) 2 -|) 

D 2 = 1 -(\q \ 2 e + 4\<lo\ 4 T 2 f+ 1 -(\q fe-l'2\<lo\ 4 r 2 ) 2 + 

(23) ^(12|g | 2 e 2 + 176Uo| 4 r 2 + l) 
Figure 1 shows the two solutions. The Peregrine and Akhmediev-Peregrine 



-IV + 
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Figure 1: (a) First-order rational solution (Peregrine breather), which ampli- 
fies the carrier-amplitude by a factor of three, (b) Second-order rational solu- 
tion (Akhmediev- Peregrine breather), which amplifies the carrier- amplitude 
by a factor of five. Both solutions are localized in both, space and time. 

breathers up to fifth-order have been recently observed in a deep-water waves 
flume p2]-[22]. We intend to conduct in the future theoretical, numerical and 
experimental investigations of NSE solutions as well as other equations of this 
kind, based on the method of simplest equation [23]- [25] and similar. 
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